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1 Introduction 



Recently, extra dimensional field theory, in particular string-derived one, plays an im- 
portant role in particle physics and cosmology When we start with extra dimensional 
field theory, it is one of most important issues how to realize a chiral spectrum in four 
dimensional (4D) effective field theory. The magnetic flux background is one of interest- 
ing ways to realize a 4D chiral theory Indeed, several field-theoretical models and string 
models, i.e. magnetized D-brane models, have been studied [UI21E1II1E1E1E1IH1I3]- 
Furthermore, magnetized D-brane models are T-duals of intersecting D-brane models. 
Within the latter type of model building, a number of interesting models have been con- 
structed m isj nni m us S 

Wavefunction profiles of zero-modes are quasi-localized on the torus with magnetic 
flux background. The number of zero-modes is determined by the size of magnetic flux. 
Since we know zero-mode profile explicitly, we can compute concretely 3-point and higher 
order couplings of 4D effective field theory by overlap integrals of zero-mode profiles [7J 
HU |15l [HI [17]. That is an important aspect of magnetized extra dimensional models. 
Moreover, such a 4D effective field theory can have Abelian and non-Abelian discrete 
flavor symmetries, which are originated from localization behavior of zero-modes in extra 
dimensions |18]J^I 

In addition to magnetic fluxes, we can introduce constant gauge backgrounds and non- 
trivial twisted boundary conditions as well as orbifold boundary conditions [7J |23j [2U [25] Jf] 
That leads richer structure in model building such as zero-mode spectra and zero-mode 
profiles. 

Non-Abelian Wilson lines, i.e. the so-called toron backgrounds [22], are also interesting 
backgrounds [301 El ED E2J EH] • They can break gauge groups with reducing their ranks. 
For a certain case with magnetic fluxes and non-Abelian Wilson lines, zero-mode profiles 
have been given [7J. Our purpose of this paper is to study more about models with 
magnetic fluxes and non-Abelian Wilson lines. We analyze zero-mode profiles in generic 
case and compute 3-point couplings. Furthermore, we study flavor symmetries. We also 
study the orbifold compactification. 

This paper is organized as follows. In section 2, we give a brief review on the extra 
dimensional models with magnetic fluxes and non-Abelian Wilson lines. In section 3, 
we study zero-mode wavefunctions on the torus compactification with magnetic fluxes 
and non-Abelian Wilson lines. We compute Yukawa couplings in section 4 and study 
flavor symmetries of our models in section 5. In section 6, we also study the orbifold 
compactification. Section 7 is devoted to conclusion and discussion. In Appendix A we 
show useful calculations, which are relevant to Yukawa couplings and in Appendix B we 
compute four-point couplings as an example of higher order couplings. 

1 See for a review |13j and references therein. 

2 Similar non-Abelian discrete flavor symmetries have been derived within the framework of heterotic 
orbifold models []j|J [501 HI] ■ Analysis on their anomalies are also important [52] . 

3 Geometrical backgrounds other than tori and toroidal orbifolds have also been studied [26j [271 HI] • 
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2 Non-Abelian Wilson lines 



2.1 Higher dimensional super Yang-Mills theory 

Our starting point is iV = 1 U (N) super Yang- Mills theory in D = 4 + 2n dimensions 
with n = 1 or 3. Its Lagrangian is written by 

£ = -^Tr {F MN Fmn) + ^Tr (AT^A) , (1) 

where M, N — 0, • • • , (D — 1). Here, A denotes gaugino fields, T M is the gamma matrix 
for D dimensions and the covariant derivative D M is given as 

D M \ = d M \-t[A M ,\], (2) 

where Am is the vector field. 

Here, we consider the torus (T 2 ) n as 2n-dimensional extra dimensions and denote their 
coordinates by y m (m — 4, • • • , 2n + 3). We use the orthogonal coordinates and choose 
the torus metric such that y m is identified by y m + 1 on the torus. The gaugino fields A 
and the vector fields A m corresponding to the compact directions are decomposed as 

H X iV) = ^Xn(x) ®Vn(y), (3) 
n 

A m (x,y) = ^2(p n>m (x) <g> 0n, m (y), (4) 

n 

where x denotes the coordinates of four-dimensional uncompact space R 3,1 . Here, we 
are interested only in zero- modes, ipo{y) and (po,m{y)- Thus, we omit the mode index 
corresponding to n = and write them as ip(y) and 4> m {y). 



2.2 Non-Abelian Wilson lines 

Here, we consider T 2 of (T 2 ) n , whose coordinates are denoted as (1/4,1/5). As a U(N) 
gauge background, we introduce the following form of (Abelian) magnetic flux, 

F„ = 2, ( '-J* l ), (5) 

where ljv a denotes (N a x N a ) identity matrix. For example, we use the following gauge, 

A A = -F452/5, A 5 = 0, (6) 
for the U(N a ) part. Then, their boundary conditions can be written as 

A m (y 4 + l,y 5 ) = A m (y^)+(^ mX ^ °) , X 4 = 0, 

A m (y 4 ,y 5 + 1) = A m (y 4 ,y 5 )+^ mX 51 ^ X5 = -27r/ a2/4 . (7) 



This background breaks the gauge group U(N) to U(N a ) x U(N — N a ). The zero mode 
ip(y) corresponding to the gaugino is also decomposed as 



A B 
C D 



depending on their U(N a ) x U (N — N a ) charges. That is, A and D correspond to the 
gaugino fields of unbroken symmetries, U(N a ) and U(N — N a ), respectively, while B 
and C correspond to bi-fundamental representations, (N a ,N — N a ) and (N a , N — N a ), 
respectively. 

Since only the U(l) part of U(N a ) has the non-trivial background, its charge q is 
relevant, that is, A, B, C and D have charges q = 0, 1, —1 and 0, respectively. For example, 
the zero-mode of B elements satisfies the following equation, 

f m (d m - iA m )B(y) = 0, (9) 

for m = 4, 5, where A m denotes the U(N a ) gauge background ([HD- Also, the zero-mode of 
C elements satisfies Y m (d m + iA m )C(y) = 0, while the zero-modes of A and D elements 
satisfy T m d m A(y) = and T m d m D(y) = 0. Here, T m corresponds to the gamma matrix 
for the two-dimensional torus T 2 , e.g. 



1 f 5 = ( -i 

10 7' ! 



(10) 



and if)(y) is the two component spinor, 



*-(£)< < u > 

that is, A, B, C and D also have two components, A±, B±, C± and D±. 

In particular, we are interested in matter fields. Thus, let us concentrate on B and C 
fields. Because of ([7]), the spinor field, e.g. B, satisfies the following boundary conditions, 

B(y A + l jy5 ) = e^B(y A ,y 5 ), (12) 
£(2/4,2/5 + 1) = e>x*B(y 4 ,y B ). (13) 

Here, we write these boundary conditions as 

5(1/4 + 1,2/5) = ^4(2/4,2/5)5(1/4,2/5), (14) 
5(2/4,2/5 + 1) = ^5(2/4,2/5)5(2/4,2/5). (15) 

The above case corresponds to ^4(2/4,2/5) = e' X4 and ^5(2/4, 2/5) — e * X5 - Then, the consis- 
tency for the contractible loop, i.e. (2/4,2/5) -> (2/4 + 1,2/5) -> (2/4 + 1,2/5 + 1) -> (2/4,2/5 + 1) -> 
(2/4,2/5) requires 

(^5 ^2/4, 2/5 + l)^ 1 (2/4 + 1, 2/5 + 1)^5(2/4 + 1, 2/5)^4(2/4, 2/5)) 5(?/4, 2/5) = 5(2/4, 2/5)- (16) 
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The left hand side reduces to e~ 27Tl ^ a ip(y 4 ,y5) in the above background. This condition 
for B leads to the quantization condition of the magnetic flux f a . That is, the magnetic 
flux f a should be quantized such that f a = integer. The consistency condition for the C 
fields also leads to the same quantization condition, i.e. f a = integer. 

When we introduce a non-trivial background for the SU(N a ) part of U(N a ), the situ- 
ation changes. That modifies the boundary conditions on, for example, B, 

B(y4 + 1, y 5 ) = y 5 )B( yi , y 5 ) = e lX4 u 4 B(y 4 , y 5 ), (17) 

B(v4,ys + 1) = ttsiVA, Vs)B(V4, Vs) = e lX5 u 5 B(y 4 ,y 5 ), (18) 

where u m are constant elements of SU(N a ). Then, the consistency condition <HB\i reduces 
to 

u}^ 1 u s u^- 2Ki f' = l Na . (19) 

If u 4 and CJ5 commute each other, that would require again e~ 2m ^ a = 1. Thus, it is 
interesting that u) 4 and CU5 do not commute each other, that is, non-Abelian Wilson lines. 
In particular, we consider the case that u) 4 oj^U} 4 corresponds to the center of SU(N a ), 
that is, 

u 5 bj A u 5 u 4 — e ijv a , {ZUj 

where M a is an integer. In this case, the consistency condition (|19p requires that the 
magnetic flux should satisfy /„ = M a /N a (mod 1). 

We denote P a = g.c.d.(M a , N a ), m a = M a /P a and n a = N a /P a E A solution of Eq. (jUJ) 
is given as 

U4 = Pa, UJ 5 = Q- m \ (21) 

where 





( 


Ip, 





o\ 




(ip. 








\ 


Pa = 








Ip. 





, Q a = P (na - 1)/2 















w 








V 




V 











(22) 



with p = e 2m l n °-. 

These non-Abelian Wilson lines break the gauge group U(N a ) further. The following 
condition on the U(N a ) gauge field, 

= w 4 A^u 4 l = w^A^ 1 , (23) 

is required. Then, the gauge group U(N a ) breaks to U(P a ). 
4 Here, g.c.d. denotes the greatest common divisor. 
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3 Matter fields 

Here, we consider the following form of U(N) magnetic fluxes, 

/ Ai Nl o \ 



'45 



2tt 



V 



o 



(24) 



/nljV„ / 



This form of magnetic fluxes breaks U(N) to Y\i U(Ni) for fa = integer. Furthermore, the 
gauge group is broken to Yli U(Pi) when we choose fa = Mi/Ni with P { = g.c.d.(Mj, iVj) 
and non-Abelian Wilson lines such that they satisfy the consistency condition like Eq. f[T5"j) . 

Now, let us focus on the (iV a + N h ) x (N a + AT b ) block in [/ (AT), which has the magnetic 
flux, 



F = 2vr 



2s i 



V„ 



We use the same gauge as Eq. 

A 4 = -2tt 



i.e. 



m b -i 



m b 2 

rif, 



2/5, 



A", 



A, = 0. 



Similarly to Eq. (j7j), we denote their boundary conditions as 

An(y4 + i,ys) = A m (y 4 ,y 5 )+( drnX l lNa R JL 
^(y*,WB + i) = An(y 4 ,y 5 )+P mX f ^ ° 6 , 

where 



(25) 



(26) 



(27) 



n a n b 



(28) 



We decompose the gaugino fields of this block in a way similar to Eq. (JSJ). That is, 
A and .D correspond to adjoint matter fields of U(N a ) and U(N b ), respectively, while B 
and C correspond to bi-fundamental representations, (N a , Nb) and (N a ,N b ), respectively. 
Among them, we concentrate on the field B, which satisfies the boundary conditions, 

B(y 4 + l,y 5 ) = ^5(y 4 ,y 5 )(^) t = e i W-^S(i/4,j/5)(^) t , 

B(y 4 ,y 5 + 1) = ngB(y4, W )(ng)t = e < W-^^S( W 4 ) y 5 )(w|) t . (29) 

Here, ui^'s are non-Abelian Wilson lines, which are given as Eqs. (l2Tj) and (|22|) . Then, 
the gauge symmetries are broken to U(P a ) and U(P b ). We study zero-mode profiles of B 
fields in what follows. 
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3.1 Integer magnetic fluxes 



Before considering the models with fractional magnetic fluxes and non-Abelian Wilson 
lines, it would be convenient to review briefly the models with integer magnetic fluxes 



and no Wilson lines, that is, 



n„ 



a,b 



a,b 



1. 



(30) 



Then, the boundary condition (|29|) reduces to 



Bp q {yi + 1,2/5) 
B pq (yi,y 5 + 1) 



B p , q (y<i,y$), 
e- 2 ^B m (y„y 5 ), 



(31) 



where m = m a — mb. Suppose that m > 0. Then, the B + component for each B v q element 
has m independent solutions for the zero-mode Dirac equation (J9]) with above boundary 
condition (131"]) . These solutions are given by 



& (2/4, 2/5) 



Z^ 1 

1 



-mw(l+^) 2 +2TTirn(l+^)y4-TTrnyl-2irm(l+^)y s 



[mz, mr), 



where z = y^+iy^, j = 0, 1, ■ ■ • , m— 1 and r — i. Here, $ 



(32) 



(mz, mr) denotes the Jacobi 



theta function. On the other hand, the B- component has no normalizable zero-modes. 
Similarly, the C_ fields have the same solutions as B + , but the C + has no normalizable 
zero-modes. 

When m < 0, the B_ and C + fields have the \m\ independent solutions with the same 
wavefunctions as above except replacing m by \m\. However, the B + and C_ fields have 
no normalizable zero-modes. 



3.2 Fractional magnetic fluxes 

Here, we study zero-mode profiles in the models with fractional magnetic fluxes and non- 
Abelian Wilson lines. 



3.2.1 



n„ 



rib 



First, let us study the magnetic flux (1231) for n = n a = n^. In this case, the non-Abelian 
Wilson lines break the gauge group U(N a ) x U(N b ) to U(P a ) x Z7(P&), where P a = N a /n 
and Pb = Nb/n. Following this breaking pattern, we decompose the fields B as 



B 



( B 00 

Bio 



\B n . 



1.0 



Bqi 

B n 

B n -i, 



\ 



(33) 
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Each of B pq components is (P a xPb) matrix- valued fields, which correspond to bi-fundamental 
(P a , Pb) fields under U(P a ) x U(Pf,). The boundary condition fl29l) due to the non-Abelian 
Wilson lines is written as 

B pq (y 4 + l,y 5 ) = B p+1>q+1 (y 4 ,y 5 ), 

B pq (y 4 ,y B + l) = p- (map - mbq) e- 2 -^y*B M (y 4 ,y 5 ), (34) 

where m is used as m = m a — mb, and B p+n>q = B Piq+n = B p _ q . That leads to the boundary 
condition, 

B pq {y 4 + n, y 5 ) = B pq (y 4 , y 5 ), 

B pq (y,,y 5 + n) = e~ 2 ™ m ^ B pq { yi , y 5 ). (35) 

Suppose that mn > 0. Then, similarly to section 13. f I the B + component for B PA has 
nm independent solutions for the zero-mode Dirac equation ([9]) with the above condition 
f l35l) . These solutions are given by 

l 

(mz,nmr), (36) 

where j = 0, 1, • • • ,nm — 1 and r = i. On the other hand, the B_ component has no 
normalizable zero-modes. One finds that these solutions satisfy the boundary conditions, 

&{y* + 1,2/5) = ^ 0^(2/4,2/5), 

e>(y A ,y B + l) = e-^& +m (y 4j y 5 ). (37) 

Thus, the zero-mode solutions with the boundary conditions (1511) due to non-Abelian 
Wilson lines can be written in terms of O- 7 of Eq. (1361) as 

n-l 

B^iyA, 2/5) = 4 q e 2m ^- m ^Q> +mr , (38) 

r=0 

where j = 0, 1, m — 1. Here, c pq is a constant normalization, which can be determined 
by the boundary conditions. Note that the boundary condition (j34|) relates 5 Pi9 and 
B p+ i^ q+ i. Thus, there are the series 

Bp,q B p+ i iq+ i B p+n q+n = B p q . (39) 

The periodicity of the series is equal to n, and there are n independent series. Each of the 
series has m independent solutions fl38|) and the total number of independent zero-modes 
is equal to mn. 

We have concentrated on the B + fields. Similarly, when mn > 0, the C_ fields have 
the same solutions as B + . However, the 5_ and C + have no normalizable zero- modes 



_vrm 2 



3 
nm 



8 



for mn > 0. On the other hand, when mn < the _B_ and C + have normalizable zero- 
modes with the same wavefunctions as the above, while B + and C_ have normalizable 
zero-modes. 

We have considered the zero- mode profiles of fermionic fields. If 4D N=l supersym- 
metry is preserved, the scalar mode has the same zero-mode profiles as its fermionic 
superpartner. 



3.2.2 n a j£ n b 

Next, we study the model with n a ^ n b . In this case, the non-Abelian Wilson lines break 
the gauge group U(N a ) x U(N b ) to U(P a ) x U(P b ), where P a = N a /n a and P b = N b /n b . 
Similarly to the previous subsection, we decompose the fields B as 



B 



( B 00 
Bio 



P>oi 
B n 



B { 



0,n b 



-1 \ 



\-Bn a -l,0 Pn a -\,\ 



(40) 



B na -l,n b -lJ 



Each of B pq components is (P a x P b ) matrix-valued fields. The boundary condition f l29|) 
due to the non-Abelian Wilson lines is written as 



BpqiVA + 1,2/5) 
Bpqiyi, y 5 + 1) 



B n 



Vf 1,<2+1 (2/4, 2/5) 



X llr, Til, ! U s~> nn P 



° q) B P , q (y. 



'4, 2/5 J 



(41) 



where B 



P+n a ,q 



B 



B p q . This boundary condition relates B p q and B p+ i ;q+ i. Then, 



similarly to f l3"§j) . there are the following series 



S M . (42) 

Here, the periodicity of the series is obtained by Q ab = l.c.m.(n , n b )$ and the number 
of independent series is equal to k ab = g.c.d.(n a ,n b ). Obviously, there is the relation, 

Q 



ab 



lib 



The above boundary condition ( 141]) leads to the boundary condition, 



Bpq{yA + QaJb,Vs) 
Bpq{yA,Vb + Qab) 



B pq (y4,y 5 ), 



2ivi j 



B p ,q{y4, 2/5) 



(43) 



Here we have defined 'the intersection number' I ab = n b m a — n a m b analogous to inter- 



secting brane models. There are S, 



ab 



n,,ni, 



I ab independent zero-mode solutions, which 



satisfy the boundary condition (1431) . Those functions are obtained as 



3 °ab ^ab °ab Q* h ^ ab °ab 



_ 7rS ab y 2 



S a b 





{{S ab /Qab)z, S ab r) 



(44) 



Here, l.c.m. denotes the least common multiple. 



9 



where t = i. These wavefunctions satisfy the following boundary conditions, 

Q j {y 4 + l,y 5 ) = e 2m ^ 3 e j ( yi ,y 5 ), 

©%4,?/5 + l) = e m( ^-^ )yi Q ] -^{y^y,). (45) 



Thus, the zero-mode wavefunctions, which satisfy the boundary conditions ( T4T]) . are ob- 
tained as 

*i(«*.wO = 4, fe^-^e^'^,^), (46) 

r=0 

where ?' = 0, 1, • • • , M — 1. Hence, the number of the independent zero-modes is equal 
to M ab = — 1^ in each of the series ( B2l . and there are the fc a b independent series. 

^cab &ab 

Thus, the total number of zero-modes is equal to M ab k ab = I a b- 

As an illustrating example, we consider the model with n a = 2, rib = 4 and m a = 
m b = 3. Then, we have k ab = g.c.d.(n a ,n b ) = 2^1, Q ab = 4, S ab = 12 and I ah = 6. We 
decompose the bi-fundamental fields B with the 2x4 matrix entries as 

B _ f Boo B 01 B 02 B 03 \ 
\Bio B u B 12 B 13 J 

From the wave function formula in Eq. (l4"6l) . one obtains the three independent solutions 
labeled by j = 0,1,2 for each component of B pq and these are represented by linear 
combinations of B l in Eq. (I44p . For example, the B 00 and B i fields are 

B 3 m = e j + e j+3 + e i+6 + e j+9 , 

B 3 01 = e J + e-^e J+3 + e- 3m e J+6 + e-^e J+9 . (48) 

Obviously, the ^-direction boundary condition can connect some of B VA components 
follows 

Eoo ->• Bu ->• 5 02 B 13 Boo, (49) 

S i ->• 5 12 5 3 B 10 -> Soi- (50) 

Since there are fc a h = 2 independent series, there are J a b = 6 zero mode solutions in this 
background. 



3.3 Another representation of solutions 

In the previous section, we have presented solutions in terms of the G- 7 functions. However, 
by using the properties of the theta function, one can represent the wave functions fl38l) 
and fj46l) as a single theta function as 

Bi pq (y 4 ,y 5 ) = C^e-^x$ 



( IabZ + i—p - —q) JabT) , (51) 
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where I a b = Iab/n a rib- The constant C 3 pq can be determined by the boundary conditions. 
The net number of zero- mode multiplicity for each of the series is given by M a b = I a b/k a b- 
Therefore the wave functions B 3 pq {y^ y 5 ) with f = j + M ab should be equal to B 3 pq {y^ y 5 ). 

p, q +n b = B 3 pq and we have twist boundary con- 
Then these conditions imply the following con- 



Furthermore we impose the B p+riaq 



dition B j (y A + 1,1/5) 



^+1,9+1(2/4,2/5) 



straint for the coefficients of C 3 as 



27ry - 



'abC J j 



°p+l,9+l 



-2m j - 



— b c 3 , =c 3 

^p,q+n b ^pqi 

(~<j+M ab _ (~tj 
p,q p,q' 



(52) 
(53) 



We start with a certain element, e.g. Cq , Then we fix other elements by using the above 
relations. If there are still unrelated elements, we start with one of those elements and 
repeat the procedure recursively again. 

In general, their solutions for C 3 pq can not be determined uniquely. Only in specific 
cases, we can write C 3 pq by a simple form. For example, when j^, = integer, C 3 pq 
is reduced to C 3 pq = const., i.e. independent of p, q and i. Most of the following models 
correspond to this case. As another example for a simple form, if we can find a certain 
integer L, which satisfies 



Mnhln 



n„ 



M ab l b + m b 



(54) 



where l a and lb are also integers, then we can write C 3 pq by the following form 



Cp q - < 



1X1 ab 



(55) 



Then the forms of wave functions would become simple as 



B 3 pq (y4,y5) 



x t) 



Mai 





hbZ + 



n„ 



-P 



ri h 



q , I ab r , (56) 



up to a normalization factor. However this expression is only valid if there exists such an 
integer L satisfying the relations. At any rate, in generic case we determine C 3 pq recursively 
as mentioned above. 

So far, we have considered the simple T 2 , where y 4 and y 5 are identified as y$ ~ y 4 + 1 
and 2/5 ~ 2/5 + I. Similarly, we can study the torus compactfication with an arbitrary value 
of the complex structure modulus r, although we have fixed r = i in the above analysis. 
Then, we obtain zero-mode wavefunctions similar to Eq. (15T|) for an arbitrary value of r 
as Eq. (15T|) . We also replace z = y i + iy 5 in the theta function by z = y^ + ry 5 . 
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4 Yukawa couplings 



Here, we study Yukawa couplings. Let us consider the following form of the magnetic 
fluxes, 



and non-Abelian Wilson lines similar to ( 12 ip . Then, there are three types of matter fields, 
(N a ,N b ), (N b ,N c ), (N c , N a ) and their conjugates under U(N a ) x U(N b ) x U(N C ), although 
they break to U (P a ) X U (P b ) x U (P c ) by non-Abelian Wilson lines. We consider the case 
with ^-^>0, ^-^>0and^-^>0. Then, the three types of matter fields 

n a n b ' n b n c n a n c ' J 1 

whose wavefunctions are denoted by tyi> Ml } \]/ fc ' M 2 ar fo ^fi,M 3 y^ a pp ear j n fa e following 
off-diagonal elements, 



where Mi = M ab} M 2 = M bc and M 3 = M ac for simplicity. We use the same indices 
for Q ab and others, i.e. Qi = Q ab , Q 2 = Qbc and Q3 = Q ac . As already explained, in 
the background with fractional fluxes and non-Abelian Wilson lines, their fields are the 
matrix valued wave functions. The Yukawa coupling can be calculated by computing the 
following overlap integral of zero-modes in the (7/4, 7/5) compact space, 



That is, the Yukawa coupling Y^ k in 4D effective theory is obtained as their products on 
(T 2 ) n , i.e. Y^ k = gn Yld=i Vd k -> wnere V % £ k denotes the overlap integral similar to Eq. (j5"9l 
for the d-th torus (T 2 ) and go is the D-dimensional gauge coupling. From this structure, 
one can see that the allowed couplings are restricted. In order to see it, we introduce 
the following parameters as k\ = g.c.d.(n a ,n b ), k 2 = g.c.d.(n b ,n c ), k 3 = g.c.d.(n a , n c ) 
and K = g.c.d.(fci, k 2 , k^) = g.c.d.(n a ,n b ,n c ). Then the parameter K determines the 
allowed couplings of Yukawa interactions. If K — 1, all of possible combinations (p,q,r) 
appear in Eq. fl59|) . However, if K ^ 1, only restricted combinations of (p,q,r) appear 
in Eq. fl59|) . but not all combinations. That is, the couplings are restricted by the 
symmetry. Indeed, allowed combinations of (p, q, r) are controlled by the gauge invariance 
before the gauge symmetry breaking. This symmetry is unbroken symmetry in the 
original gauge symmetry. 

Now, let us consider the following summation of wavefunction products, 




(57) 




(58) 




(59) 



rjkl 
pgr 



•Sji \If fc filr' )* 4- vT/J 

pq qr\ rp/ ' p 
+ • • • + \I/„ 1 _i ,r 




(60) 
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where Q = l.c.m.(Qi, Q2, Qz)- One can represent Q as Q = Q\q\ = Q 2 q2 = QsQs- To 
compute the integral it is useful to represent the wavefunctions as follows 





k' ' 



' v_ 

M! A 





M[ fm a m b \ M[ 

—z + — p q , —r , 

Q \n a n b J Q 



M' 2 



Q 



m b m c 



—z+ — q r ,—t , 



Q 



n b n c 



m a m 



M' 2 



— >-z + — p-— r ,— >-t , 



n r 



Q 



M> 



(61) 



Q 



where j' = qxj, k' = q 2 k, I' = q 3 l and M[ = qiMi, (i = 1,2,3). Here the relation 
M[ + M 2 = M3 holds. By using the product property of the theta function, the product 
of \$fjMi\$jk,M2 j s represented by the sum of the theta functions as 



PI 1 r pq qr 



+k'+M[m 



meZ, 



M3 

~Q 



■z + 



in,, m c \ \ /' 
P r 5 TT T 



n„ 



Q 



xtf 



M^j'-M[k'+M[M^m 
M[M' 2 M' 3 



— M'p -M'q -M[q + —M[r, 1 2 3 t 

n a n b n 6 n c Q 



(62) 



Here one can use the properties of boundary conditions for non-Abelian Wilson lines. 
Using the property of V p , q (y4 + l,y 5 ) = ^p+1,,+1 (2/4,2/5) we find 



/•l rl i-Q rl 

/ dy A / dy 5 r£ = dy A dy 5 % q ¥ qr (V 
Jo Jo Jo Jo 



k r P r- (63) 



Therefore we can obtain the analytic form of Yukawa couplings and flavor structures 
similar to the case with Abelian Wilson lines. By using the orthogonal condition for the 
matrix valued wave functions (see Appendix A) as 



f Q dy 4 f ' dy^^f^ = 6 jlk , 
Jo Jo 



(64) 



one can lead to the following form of Yukawa couplings 



£ dy, jf 1 dy 5 % q V qr (* k rp r = Q J E 



Sj'+k'+M^m, J'(modM£) 



X d 



M^j'~M[k'+M[Mf 2 m 

WWW 




m a ~ 



m b 



Q — hcP H h a q H hbr , 



n„ 



n b 



M[M' 2 M' Z 

~Q 



.(65) 



up to the factor N Ml N M2 N* M3 Ci q Cl r (C l pr y . Here, the Kronecker delta 5 r+kl+M ^ ;/( mo dM^ 
leads to the coupling selection rule 



f + k' + M[m = I' mod M' 3 , 



(66) 
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where m = 0,1,...,M' 3 - 1. When g = g.c.d.(M(, M' 2) M3) = g.c.d.(Mi, M 2 , M 3 ), the 
coupling selection rule is given by 

j' + k' — I' mod (67) 

That means that we can assign Z g charges to all of zero-modes]^ 

Here we study again the Z% symmetry, which we showed. The total number of multi- 
plicity of is nothing but |I a &|, and it is represented by two parameters of k a b and M a b 
as I a b = k a bM a b. If K = g.c.d. (k a b, kb c , k ca ) 7^ 1, they are divided to K types of zero-modes 
and distinguished by labeling the component of each matrix. We introduce such a kind 
of flavor indices as j, k and / for ab-, be-, ca-sectors, respectively. We define the relation 
between the flavor labeled by j and the component of matrix p,q as j = p — q mod k\. 
Similarly the other sectors are also defined as k = q — r mod &2 and I = p — r mod k 3 . 
Since the allowed couplings must be gauge invariant, there is the coupling selection rule 
for this kind of flavor indices, which is given by 

j + k = I mod K. (68) 

This is because the Yukawa couplings are restricted in the trace of the matrix. Therefore 
we find two types of coupling selection rules, i.e. the Z g and Z K symmetries. 

We can extend the computation of 3-point couplings to higher order couplings. For 
example, in appendix B, we show the computation of 4-point couplings. 



5 Non-Abelian discrete flavor symmetry 

Here, we study the non-Abelian flavor symmetries, which can appear in our models. 
5.1 The case with M{ 7^ 1 and k{ = 1 

First, we consider the models with k\ = k 2 = k 3 = 1. Then, the number of zero- 
modes are given by \I a b\ = M±, |/& c | = M 2 and |J ca | = M 3 . We consider the models 
with g = g.c.d.(Mi, M2, M3) ^ 1. The Yukawa couplings do not depend on the matrix 
components {p,q,r), and are reduced to the following form 



j'+k'+M[m,l'(modM' 3 ) 

1 

3 

M' 2 f '~M[k '+M[ML,m ' 







(0,M(M^M^/Qr), (69) 



where we have taken simply p = q = r = and the phase factor like C 3 pq disappears. This 
form is nothing but the case with integer fluxes and without non-Abelian Wilson lines. 



6 See Refs. [331 133 for a similar selection rule in intersecting D-brane models. 
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In this types of Yukawa couplings, 4D effective theory has another flavor symmetry called 
by the shift symmetry, which corresponds to the transformations of flavor indices as 



j'^j' + M[/g, 
k '^k' + M' 2 /g, 

+ M^/g, 



(70) 



simultaneously. Under this transformation, Yukawa couplings are invariant. There is 
also coupling selection rule as shown in the previous section given by the Z g symmetry 
(1671) . Then, they form the non-Abelian discrete flavor symmetries as the same as the case 
without non-Abelian Wilson lines. 

For simplicity, suppose that M[ = g. Then, there are g zero-modes of \E r3 ' ,Ml . The se- 
lection rule ( 1671) means that 4D effective theory is symmetric under the Z g transformation, 
which acts on tyi ,9 as Z^i ' s , where 



/ I 



\ 



V 



(71) 



P 9 - 1 ) 



and p = e lm l g ' . Furthermore, the effective theory has another symmetry flTOj) . That can 
be written as cyclic permutations on ,g , 



3 ,9 



(72) 



That is nothing but a change of ordering and also has a geometrical meaning as a discrete 
shift of the origin, z = — > z = — -. This symmetry also generates another Z g symmetry, 

(C) 

which we denote by Z g and its generator is represented as 



C 



( 1 
10 



o\ 





(73) 



\ 1 ••• J 

on \f/i'' 9 . These generators, Z and C, do not commute each other, i.e., 

CZ = pZC. 



(74) 



Then, the flavor symmetry corresponds to the closed algebra including Z and C . Diagonal 
matrices in this closed algebra are written as Z n (Z') m , where Z' is the generator of another 



Z' g written as 



Z' 



(p 

V 



p) 



(75) 
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w 


Representation of D4 


2 


2 


4 




6 


3x2 



Table 1: _D 4 representations of zero-modes in the model with g = 2. 





Representation of A (27) 


3 


3 


6 


2x3 


9 


li, I2, I3, I4, I5, 3-6; 3-7, 3. 8 , lg 


12 


4x3 


15 


5x3 


18 


2 X 1 2 , I3, I4, I5, 16, I7, 3.8; 3-9} 



Table 2: A(27) representations of zero-modes in the model with g = 3. 

on vl/- 7 ' 9 . Hence, these would generate the non-Abelian flavor symmetry \Z g x Z g ) x Z g , 
since Z 9 x is a normal subgroup. These discrete flavor groups would include g 3 elements 
totally. 

For example, for g = 2 and 3 these flavor symmetries are given as ^2x1^2 = -D4 and 
(Z3 x Z3) x Z3 = A(27), respectively. Then, the fields ' 9 correspond to 2 of -D4 and 3 
of A(27), as shown in Tables [1] and [2], respectively. When M'/g is an integer larger than 
1, the ^ ' fields correspond to other representations. For smaller values of M'/g, the 
corresponding representations are shown in Tables [T] and [2J 

However we note that their multiplets have several types of representation under this 
symmetry. Because a Z g charge of fields labeled by j is not j but j' = qj. Therefore even 
if they have same multiplicities (Mi = M2), their representations may be different from 
each other. 

5.2 The case with M % = \ and k ^ 1 

Next, we consider the models with = 1 and k ^ 1. In this case, we also find flavor struc- 
tures similar to the case without non-Abelian Wilson lines. Suppose all the components 
of zero modes are given by \I a i\ = hi, \Ib c \ = k 2 and |J co | = k 3 . Then it is possible to take 
phase factors for each of wave functions C 3 pq = 1. We commonly use K = g.c.d.(fci, k 2 , k 3 ). 
The Yukawa couplings depend on the indices p, q and r only through a combination 9 pqr 
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given by 



a n m «f , mb f , mc f 

\n a n b n c 

= Q ( —hc(j + nikx) - —I a b(k + n 2 k 2 ) J , (76) 

where we have used the relations p — q = n\k\ + j and I — r — n 2 k 2 + with ni,n 2 G 
We find that the Yukawa couplings are invariant under the following transformation as 



3^3 + 



K 



k + (77) 



It is obvious that this transformation is the permutation of flavor index with order K. 
Therefore we have two symmetries: one is the discrete Zk symmetry comes from the 
coupling selection rule and another is this shift symmetry. By combining these two sym- 
metries, it become the same non-Abelian discrete flavor symmetry as the case without 
Non-Abelian Wilson-lines. That is, these flavor symmetries are given as Z 2 x Z 2 = _D 4 for 
K = 2, (Z 3 x Z 3 ) x Z 3 = A(27) for K = 3 and (Z K x Z K ) x Z K for generic K. 

We have two aspects of flavor structures which are characterized by the parameters 
M, K. In the latter case, the origin of flavor symmetry is the gauge symmetry. The 
background breaks the continuous gauge symmetry, but discrete symmetry remains as 
the flavor symmetry. In the former case, the flavor would not directly originated from the 
gauge symmetry. However, T-duals of both cases would correspond to similar intersecting 
D-brane models, where n a and m a have almost the same meaning, that is, winding num- 
bers of D-branes for different directions. Thus, these two pictures of flavor symmetries 
are related with each other by T-duality through the intersecting D-brane picture. 

So far, we have considered the models with Mj = 1 and K ^ 1 and found the flavor 
symmetry (Zk x Zk) x Z k . Here we comment on generic case with M ^ 1 and K ^ 
1. Even in such a case, the selection rules due to Z g and Z^ symmetries hold exact. 
However, the general formula of Yukawa couplings depend on both the indices j and j. 
Then, 4D effective Lagrangian is not always invariant under the above (independent) shift 
transformations (!70|) and (J77J). 



5.3 Illustrating examples 

We show two illustrating examples. We concentrate on only the T 2 torus. The first 
example is the model with (Ii,I 2 ,I 3 ) = (2,4,2). The background magnetic flux is taken 

as 
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Then the appearing chiral matters are denoted by 



const Li' Ml_1 



A 



ffl,M 3 =l 



rp 



PI 

const R^ M2=1 | , 
const 



(79) 



where p =0, 1, q =0, 1, 7 and r =0, 1, 2, 3. The wave functions are represented by 
following theta functions as 



LUy) 



R k q M 



HUy) = Nu^l^-d 



(z/S + (l/2p-3/8q),r/S), 
(z/8 + (3/8q-l/Ar),r/8), 
(z/4 + (l/2p-l/4r),2r/8), 



(80) 



where we take j = k = I = 0. The several parameters are also given by these fluxes. We 
have k\ = 2, k 2 — 4, k 3 = 2 and K = g.c.d.(/ci, k 2 , k 3 ) = 2. The gauge invariant 3-point 
couplings are divided to four types of Yukawa couplings shown below 



C 

£-000 

Am 
£101 



^000 + ^OlO + A)01 + £l0(b 

LqqRooHqq + LiiRnH n + Lq 2 R 22 H 02 + -^13-^33-^13 
+L 04: R4qH q + L 15 i? 51 i^ 11 + L 06 Rq 2 H 02 + Li 7 R 73 H 13 , 
LooRoiH 01 + L\\R\ 2 H 12 + Lq 2 R 23 H 03 + Li 3 R 3 qH w 

+-^04-^41-^01 + -^15-^52-^12 + -^06-^63-^03 + -^17-^70-^10? 
L\qRqoH\q + Loi-Rll-f^Q! + L\ 2 R 22j H\ 2 + -^03-^33-^03 
+Li4i?40-ff| + Lq^R^iH^ + LiqRq 2 H\ 2 + LqjRjsH^, 

LwRoiH^ + L i-Ri2-f^oi + Li 2 R 23 H\ 3 + L 03 R 30 Hq 
+L u R^iH\ 1 + L 05 R 52 Hq 2 + L\ & R &3 H\ 3 + LqyR 7 qHq Q . 



As seen in these interaction terms, one finds that all the combinations (p, g, r) are not 
allowed. This is because it has = g.c.d.(2, 4, 2) = 2. These fields L,R,H are divided 
to two classes under the discrete Z 2 charge. For instance, for R fields, the flavor index is 
defined by k = q — r mod 4. We assign the Z 2 charges as 



fc=0 D fc=2 



Z 2 + : iT =u , R 



)fc=l D fc=3 



Zj 2 — .It , ix 

and for other fields we also assign the Z 2 charges as 



(81) 



Z 2 + 
Z 2 - 



£fc=0 jjk=0 
L k=l H ~k=l_ 



(82) 
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That corresponds to the coupling selection rule as j+k = I mod 2. The Yukawa couplings 
Yj™ are obtained after the overlap integrals as 



Yff oc d 



(l/2p - 3/4g + l/4r, r/4). 



(83) 



We also consider about the shift symmetry for this model, i.e. 

~ m c I ab 



3^3 



K 



— j + 1 mod 2, 



k -> fc + !!^ = fc + 2 mod 4, 

[ [ + H^ = [ + i mod2 . 



(84) 



As shown in the previous section, the Yukawa couplings are invariant under this transfor- 
mation. These two operators make the D A = Z 2 x Z 2 discrete flavor symmetry. One can 
understand the representation for each field under D 4 symmetry. As an analysis similar 
to the previous section, one can find that L and R correspond to doublets and H fields 
become four non-trivial singlets under Z) 4 symmetry. 

As another example, we consider the model with (7 1 ,/ 2 ,7 3 ) = (3,3,3), which is not 
realized by only integer fluxes. We choose fluxes as 



'31 



N a 



F = 2tt 



2*-N b 



35) 



01, 



Then the appearing chiral matter fields are denoted as follows, 



A 



const L J 0p Ml=3 



const R 



o-«,M 3 =3 



const 



(86) 



where p, q = 0, 1. This model has fcj = 1 and Qi = 2, Q 2 = 2, Q% = 1, Q = 2 (f = j, k' = 
k,l' = 21). The gauge invariant 3-point couplings are given as 



tvLpq Rqr Hp r 



LqqRoqHIq + Lqi-RioT/qq. 



(87) 



The Yukawa couplings Yj® are calculated by overlap integrals as follows 



yjkl 
pqr 



'0 



dy 4 I dy 5 L pq (y)R k qr (y)H l rp ( y y 

(0,27r), 



o 



OC ^ 8j> + k'+2,m,l''& 



3j'-3k'+9m 
54 





(88) 
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where we take p = q = r = 0. From the structure of Kronecker delta, one can read the 
selection rule as 



f + k' + 3m = I' mod 6 
-> j + k - 21 = mod 3 



(89) 



Since is defined by g = g.c.d.(M(, M' 2 , = 3, this model has A(27) = (Z 3 x Z 3 ) x Z 3 
flavor symmetry. Here we mention that the charge assignment is different from the case 
with Abelian Wilson line. For H fields, their Z 3 charges are obtained as V = 21, and 
they correspond to the multiplet of 3 representations. Other fields, L and R, correspond 
to 3 representations, and they can couple in the language of flavor symmetry. Therefore 
the extension to the non- Abelian Wilson line case causes to have more various types of 
representations and flavor structures. 

It is possible to introduce the constant gauge potential called by the Abelian Wilson 
line. We take the previous model with (Ji, I 2 , ^3) = (3, 3, 3). We assume N a — 4, N b — 4, 
N c = 2. Then the fractional fluxes with non-Abelian Wilson lines can reduce the rank 
of gauge symmetry, that is, the £/&(4) gauge group breaks to L^(2) and the total gauge 
symmetry is U(4) a x U(2) b x U(2) c . To break the gauge symmetry U(4) a x U (2) b x U{2) c to 
the standard-model gauge group, Abelian Wilson lines can be introduced. For example, 
we can introduce the Abelian Wilson lines in U(4) a along the following direction, 



Then, the gauge group U(4) a is broken to U(3)xU (1). Similarly, we introduce the Abelian 
Wilson lines in U(2) c along the following direction, 



Then, the gauge group U{2) c is broken to U{1) x U{1). Then the (supersymmetric) stan- 
dard model with three generations is realized up to U(l) factors. We can also introduce 
the Abelian Wilson line along the U(l) direction of Z7(2)&. Since the different Wilson 
line leads to different Yukawa couplings, that would lead to various flavor structures. For 
example, the above model leads to the A(27) flavor symmetry in generic values of Wilson 
lines as studied in the previous section. However, the flavor symmetry is enhanced to 
the A (54) symmetry when Wilson lines vanish. Thus by choosing the particular choice 
of Abelian Wilson lines, we could realize that the flavor symmetry is large like A (54) in 
a subsector, e.g. in the lepton sector, but the other sector, e.g. the quark sector, has the 
smaller flavor symmetry like A(27)J3 

7 Indeed, non-Abelian discrete flavor symmetries such as D4, A (27) and A (54) would lead to phe- 
nomenologically interesting models [36l EH [38] . 




(90) 




(91) 
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6 Magnetized orbifold background 



We have obtained the explicit wavefunctions on the torus background. Here, we study 
about the models on the orbifold background. Following Ref. [23], we study the T 2 jZ 2 
orbifold, which is constructed by dividing T 2 by the Z 2 projection z — > —z. Furthermore, 
we require the field projection of periodic or anti-periodic boundary conditions consistent 
with the Z 2 orbifold, 

*H/4, -y 5 ) =P*{y4,y 5 ), (92) 

where P is +1 or — 1. One can show that the matter wave functions satisfy the following 
property 

H q (-VA, -y 5 ) = *:l- g (2/4, ?/b)- (93) 

For the case with k — 1, this relation holds, because every sector of (p,q) are related by 
the boundary conditions, so the labels (p, q) have no meaning. However, in the k 7^ 1 
case, they have k x M independent zero-modes and we symbolically denote them by ^Z- 7 '- 7 
(j = 0, 1, M — 1 and j = 0, 1, k — 1). For example, in the case with n a = rib = 3, we 
may use the following notations 

¥ m , *i, ¥ 22 -> viz'-.' °, 

^01, *L *20 -> '• (94) 

where j = p — q mod i^. Then, the above property ( |93i) can be written as 

* iJ (-2/4, -2/5) = * //.-,). (95) 

Then the even and odd wave-functions are easily obtained. For the case with M = 3, 
there are 3x3 independent fields and they are divided into the following even and odd 
wave functions 

even : tf 1 - + V 2 ' , V ' 1 + V°> 2 , V 1 ' 1 + ^ 2 ' 2 , V 2 ' 1 + V 1 ' 2 , 

odd : W 1 ' - * 2 ' , tf 1 ' 1 - ^ 2 ' 2 , q 2 > 1 - q 1 ' 2 . (96) 

Note that these represent the wave functions e.g. ^\ 2 + V&li by ^ r * L ' 1 + \l/ 2 ' 2 . As examples, 
the zero-mode numbers of even and odd wave functions for smaller values of k and M are 
shown in Table [3J 

Yukawa couplings as well as higher order couplings can be computed on the orbifold 
background by overlap integrals of wavefunctions in a way similar to the torus models. 
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k = 1 even 
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k = 2 even 
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odd 
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1 


2 2 


odd 
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k = 3 even 
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10 


k = 4 even 


3 


6 
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10 


11 


14 


odd 


1 


2 


4 


5 


7 


8 


odd 


1 


2 


5 


6 


9 


10 



Table 3: The numbers of even and odd zero-modes 

7 Conclusion 

We have studied the flavor structure of 4D effective theories, which are derived from 
extra dimensional theories with magnetic fluxes and non-Abelian Wilson lines. We have 
obtained zero-mode wavefunctions for generic case. Their Yukawa couplings as well as 
four-point couplings have been computed. Furthermore, we have also studied non-Abelian 
flavor symmetries and some parts of them are originated from gauge symmetries. In 
addition, the orbifold compactification has been discussed. 

We have obtained quite rich flavor structure compared with the magnetized torus/orbifold 
models with Abelian flavor structures. For example, there are various ways of model build- 
ing leading to three generation models. Thus, it would be interesting to apply our analysis 
to phenomenological model building. 
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A Integral of wave functions 

In this appendix we show some calculations on integral of wave functions, which are used 
in section 4. The general wave functions are expressed as 

(97) 
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n„ 



-p 



m b \ M' 
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where Q is given by integer satisfying the relation Q = k x l.c.m.(n a , n b ), k E Z. The 
orthogonal condition implies 

[ Q dy, fdy&ig (* k p > M Y = |tf*rqk(C&r /% 4 fdy^ M l^ 
Jo Jo Jo Jo 

x e _ 7r Af'/Qa+j"A/') 2 e 2 -^('+j"/A^)(?/4+ij/5+(^ P -^ 9 )) 

x ^ e --M7Q( m +i7M') 2 e -2«^(rn+j'/Af')(?/4-^5+(^p-^g)) > ^ 



Here the integral over 2/4 is obtained as 



m(mod M)$j, k(mod M) 



Then, we obtain 



(99) 



dy, I dy 5 V> q M ^ k pq M y 



.-, 1 2wM „,2 



j,fc(mod M) 



(100) 



By this we can fix the normalization. 

B N-point coupling 

One can calculate the N-point couplings by using the addition formula of the theta func- 
tions as in Ref. [15J. For example, we show here the explicit calculation for general four 
point couplings. We assume that I ab , I bc , led > and I da < 0. Four zero-mode wavefunc- 
tions are written as 
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where Q is defined as Q = l.c.m.(n a , rib, n>c, no). First, the product of ip^' q Ml and ip^ 2 
becomes 



<^; A/2 = c Pq c qr e-^ £ * 

x$ 



j'+k'+M[m 
IP 





M' 



,m a 



M' 
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where M' = M[ + M 2 . Then we repeat this product for ipl^ 3 an d use the orthogonal 
condition for the M' 4 sector because the relations M[ + M 2 + M3 = M' + M3 = hold 
by definition. Finally we obtain the overlap integral for four wave functions as 
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This result is just the product of two theta functions. By solving the Kronecker delta, we 
obtain the sum of two theta functions like ^2 rn k m y 1 m ■ Therefore even including the 
non-Abelian Wilson lines we obtain results which are similar to Ref. [TJ] for general four 
point couplings. 
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